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^ Abstract. Current status of theoretical researches on the QCD phase diagram 

QQ at finite temperature and baryon chemical potential is reviewed with special 

y—^ emphasis on the origin of various phases and their symmetry breaking patterns. 

Topics include; quark deconfinement, chiral symmetry restoration, order of 

I I the phase transitions, QCD critical point (s), colour superconductivity, various 

f | inhomogeneous states and implications from QCD-like theories. 
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1. Introduction 

One of the most crucial properties in the non-Abelian gauge theory of quarks and 
gluons, Quantum Chromodynamics (QCD) [1], is the asymptotic freedom [2, 3]; the 
coupling constant runs towards a smaller value with increasing energy scale. It is 
hence a natural anticipation that QCD matter at high energy densities undergoes a 
phase transition from a state with confined hadrons into a new state of matter with 
on-shell (real) quarks and gluons. 

There are two important external parameters for QCD in equilibrium, the 
temperature T and the baryon number density n^. (In the grand canonical ensemble, 
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the quark chemical potential /Xq — (Ib/S may be introduced as a conjugate variable 
to the quark number density n^ = Sne-) Since the intrinsic scale of QCD is 
Aqcd ~ 200 MeV, it would be conceivable that the QCD phase transition should 
take place around T ^ Aqcd ~ 0(10^^) K or tt-b ^ ^qcd ^ ^ fm^'^. Experimentally, 
the heavy-ion collisions (HIC) in laboratories provide us with a chance to create hot 
and/or denser QCD matter and elucidate its properties. In particular, the Relativistic 
Heavy-Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) has conducted 
experiments to create hot QCD matter (a quark-gluon plasma or QGP in short) by 
the An- An collisions with the highest collision energy y/s~^ — 200 GeV [4] . The Large 
Hadron Collider (LHC) at CERN will continue experiments along the same line with 
higher energies [4, 5]. Exploration of a wider range of the QCD phase diagram with 
riB up to several times of the normal nuclear matter density uq ~ 0.17 fm~ may be 
carried out by low-energy scan in HIC at RHIC as well as at the future facilities such 
as the Facility for Antiproton and Ion Research (FAIR) at GSI, the Nuclotron-based 
Ion Collider Facility (NICA) at JINR and the Japan Proton Accelerator Research 
Complex (J-PARC) at JAERI. 

In nature, the deep interior of compact stellar objects such as neutron stars would 
be the relevant place where dense QCD matter at low temperature is realized (see [6] 
for a review). In fact, continuous efforts have been and are being conducted in the 
observations of neutron stars to extract information of the equation of state of dense 
QCD. If the baryon density is asymptotically high, weak coupling QCD analyses 
indicate that the QCD ground state forms a condensation of quark Cooper pairs, 
namely the colour superconductivity (CSC). Since quarks have not only spin but also 
colour and flavour quantum numbers, the quark pairing pattern is much more intricate 
than the electron pairing in metallic superconductors. 

In this review we will discuss selected topical developments in the QCD phase 
diagram with an emphasis on the phases at finite /in. For the topics not covered in 
the present article, readers may consult other reviews [7, 8, 9, 10, 11, 12, 13, 14, 15, 
16, 17, 18, 19]. 

We organize this article as follows. In section 2, after a brief introduction of 
two key features of hot/dense QCD, i.e. the deconfinement and chiral restoration, we 
show a conjectured QCD phase diagram in ji^-T plane. In section 3 we introduce 
three major order parameters to characterize QCD matter at finite T and /xb, i.e. the 
Polyakov loop, the chiral condensate and the diquark condensate. In section 4 we 
summarize the current status on the chiral phase transition at finite T with /ub = 0. 
In section 5 we go into more details about the phase transitions in the density region 
up to a moderate /iB comparable to T, which are classified according to different 
spontaneous chiral symmetry breaking (S^B) patterns. We then address the situation 
with /xb much greater than T in section 6. Even though the theoretical understanding 
has not been fully settled down yet, section 7 is devoted to a review over conjectured 
inhomogeneous states of QCD matter as well as some alternative scenarios. In section 8 
we look over suggestive results and implications from QCD-like theories. Section 9 is 
devoted to the summary and concluding remarks. 

2. QCD phase structure 

We will run through the QCD phase transitions and associated phase structure here 
before detailed discussions in subsequent sections. 
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2.1. Deconfinement and chiral restoration 

A first prototype of the QCD phase diagram in T-ub plane was conjectured in [20]. 
It was elucidated that one could give an interpretation of the Hagedorn's limiting 
temperature in the Statistical Bootstrap Model (SBM) [21] as a critical temperature 
associated with a second-order phase transition into a new state of matter. The weakly 
interacting quark matter at large ub due to the asymptotic freedom had been also 
recognized [22]. Historical summary of QCD phase diagram and its exploration in 
HIC experiments are given in the reviews [23, 24]. 

• Deconfinement — In an early picture of hadron resonance gas at finite 
temperature [25], the density of (mostly mesonic) states p(m) as a function of the 
resonance mass m, is proportional to exp(m/r^) where T^ ~ 0.19 GeV is known 
from the Regge slope parameter. Such an exponentially growing behaviour of 
the density of states should be balanced by the Boltzmann factor exp(— m/T) in 
the partition function. When T > T^, the integration over m becomes singular, 
so that T^ plays a role of the limiting temperature (Hagedorn temperature) 
above which the hadronic description breaks down. This argument is applied 
to estimate the critical value of /ie as well. The density of baryonic states, 
Pb(wb) oc exp(?TiB/rj^), is balanced by the Boltzmann factor exp[— (tob— /iB)/r], 
leading to the limiting temperature T = (1 — fiB/mB)T^. We see equivalently 
that the critical /iB at T = is given by itib ( ^ 1 GeV) . 

If one bears in mind a simple bag-model picture that hadrons are finite-size objects 
in which valence quarks are confined, it is conceivable to imagine that hadrons 
overlap with each other and start to percolate at the Hagedorn temperature 
[26, 27], which is an intuitive portrayal of quark deconfinement. Although the 
Hagedorn/percolation picture is useful for practical objectivization, it is necessary 
to develop a field-theoretical definition of the quark-deconfinement in QCD. 
Global centre symmetry of pure gluonic sector of QCD gives such a definition 
as elucidated in section 3.1. 

• Chiral restoration — The QCD vacuum should be regarded as a medium with 
full of quantum fluctuations that are responsible for the generation of non- 
perturbative quark mass. In hot and dense energetic matter, quarks turn bare 
due to asymptotic freedom. Therefore, one may expect a phase transition 
from a state with heavy constituent quarks to another state with light current 
quarks. Such a transition is called chiral phase transition named after the 
underlying chiral symmetry of QCD. The QCD phase diagram at finite T 
and fj,B was also conjectured from the point of view of chiral symmetry [28]. 
In this case, the order parameter is the chiral condensate (iptp) which takes a 
value about — (0.24GeV)^ in the vacuum and sets a natural scale for the critical 
temperature of chiral restoration. In the chiral perturbation theory (xPT) the 
chiral condensate for two massless quark flavours at low temperature is known 
to behave as (V5V)t/(V'V') = 1 - T^/i^f^) - T*/{38Af^) - ■■■ with the pion 
decay constant /^ ~ 93 MeV [29] . Although the validity of xPT is limited to low 
temperature, this is a clear evidence of the melting of chiral condensate at flnite 
temperature. At low baryon density, likewise, the chiral condensate decreases as 

{'ipip)n^/{ip^l^) = 1 - cr^rw riB/C/^TO^ ) [30, 31, 32] where a^N - 40 MeV is the 

7r-7V sigma term. (For higher-order corrections, see [33, 34].) 
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Figure 1. Conjectured QCD phase diagram with boundaries that define various 
states of QCD matter based on SxB patterns. 



The chiral transition is a notion independent of the deconfinement transition. In 
section 3.2 we classify the chiral transition according to the S^B pattern. 



2.2. Conjectured QCD phase diagram 

Figure 1 summarizes our state-of-the-art understanding on the phase structure of QCD 
matter including conjectures which are not fully established. At present, relatively firm 
statements can be made only in limited cases - phase structure at finite T with small 
baryon density (/ib <C T) and that at asymptotically high density (/ie S> Aqcd)- 
Below we will take a closer look at figure 1 from a smaller to larger value of ns in 
order. 



Hadron-quark phase transition at ^b = 0^ The QCD phase transition at finite 
temperature with zero chemical potential has been studied extensively in the numerical 
simulation on the lattice. Results depend on the number of colours and flavours as 
expected from the analysis of effective theories on the basis of the renormalization 
group together with the universality [35, 36]. A first-order deconfinement transition 
for TVc = 3 and iVf = has been established from the finite size scaling analysis 
on the lattice [37], and the critical temperature is found to be T^ ~ 270 MeV. For 
Af > light fiavours it is appropriate to address more on the chiral phase transition. 
Recent analyses on the basis of the staggered fermion and Wilson fermion indicate a 
crossover from the hadronic phase to the quark-gluon plasma for realistic u, d and s 
quark masses [38, 39]. The pseudo-critical temperature Tpc, which characterizes the 
crossover location, is likely to be within the range 150 MeV — 200 MeV as summarized 
in section 4.2. 

Even for the temperature above Tpc the system may be strongly correlated and 
show non-perturbative phenomena such as the existence of hadronic modes or pre- 
formed hadrons in the quark-gluon plasma at /iB = [28, 40] as well as at /xb ^ 
[41, 42, 43]. Similar phenomena can be seen in other strong coupling systems such as 
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Figure 2. Characteristic points on the QCD phase diagram. E represents so- 
called the QCD critical point. F is another critical point induced by the quark- 
hadron continuity. G is the critical point associated with the liquid-gas transition 
of nuclear matter. H refers to a region which looks like an approximate triple 
point. See the text for details. 

the high-temperature superconductivity and in the BEC regime of ultracold fermionic 
atoms [44]. 

QCD critical points: In the density region beyond fiB ^ T there is no rehable 
information from the first-principle lattice QCD calculation. Investigation using 
effective models is a pragmatic alternative then. Most of the chiral models suggest that 
there is a QCD critical point located at (/ie — IJ-EtT — Te) and the chiral transition 
becomes first-order (crossover) for /ie > /iE (fJ-B < Me) for realistic u, d and s quark 
masses [45, 46, 47, 48] (see the point E in figure 2). The criticality implies enhanced 
fluctuations, so that the search for the QCD critical point is of great experimental 
interest [49, 50]. 

There is also a possibility that the first-order phase boundary ends at another 
critical point in the lower-T and higher-/iB region whose location we shall denote by 
(/^f,2f) as shown by the point F in figure 2. As discussed in section 6, the cold 
dense QCD matter with three degenerate fiavours may have no clear border between 
superfluid nuclear matter and superconducting quark matter, which is called the quark- 
hadron continuity. 

In reality, the fate of the above critical points (E and F) depend strongly on the 
relative magnitude of the strange quark mass rus and the typical values of T and /ie 
at the phase boundary. 

Liquid-gas phase transition of nuclear matter. Since the nucleon mass is mpf ~ 
939 MeV and the binding energy in isospin-symmetric nuclear matter is around 
16 MeV, a non- vanishing baryon density of nuclear matter starts arising at /xb = 
/Unm — 924 MeV at T = 0. At the threshold fis = Mnm, the density ne varies from 
zero to the normal nuclear density ng = 0.17 fm~ . For < rie < no the nuclear 
matter is fragmented into droplets with ne = uq, so that tt-b < uq is achieved on 
spatial average. This is a typical first-order phase transition of the liquid-gas type. 
The first-order transition weakens as T grows and eventually ends up with a second- 
order critical point at (/xg, Jg) as indicated by the point G in figure 2. Low energy 
HIC experiments indicate that ^q ^ ^nm and Tq = 15 ^ 20 MeV [51]. 

Quakyonic matter. The Statistical Model is successful to reproduce the experimen- 
tally observable particle abundances at various ^s„„ and thus various /xb- This 
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model assumes a thermally equilibrated gas of non-interacting mesons, baryons and 
resonances for a given T and /ie [52, 53, 54]. Within the model description one can 
extract T and /ib from the HIC data by fitting the particle ratios. The accumulation 
of extracted points makes a curve on the n^-T plane, which is called the chemical 
freeze-out line. 

The freeze-out line is not necessarily associated with any of QCD phase 
boundaries. Nevertheless, there is an argument to claim that the sudden freeze- 
out of chemical compositions should take place close to the phase transition [55]. 
Along the freeze-out line the thermal degrees of freedom are dominated by mesons for 
AtB ^ niN- The higher ^b becomes, the more baryons are excited. This indicates that 
there must be a transitional change at (TH,/iH)7 where the importance of baryons in 
thermodynamics surpasses that of mesons. This happens around /xh — 350 ^ 400 MeV 
and Th = 150 ^ 160 MeV according to the Statistical Model analysis. 

It is interesting that such a phase structure is suggested from the large N^ limit 
of QCD. When TVc is large, quark loops are suppressed by l/Nc as compared to 
gluon contributions [56, 57]. A finite baryon number density arises and the pressure 
grows of order Nc once ^-q becomes greater than the lowest baryon mass M-q . Such 
cold dense matter in the N^ = oo world is named quarkyonic matter [58]. Then, 
the phase diagram of large- iVc QCD consists of three regions separated by first-order 
phase transitions, i.e. the confined, deconfined and quarkyonic phases. The meeting 
point of the three first-order phase boundaries is the triple point whose remnant for 
finite Nc is indicated by the point H in figure 2, as suggested in [59]. We will revisit 
the idea of quarkyonic matter in section 8.1. 

Colour superconductivity: If /iB is asymptotically large, i.e. /xb S> Aqcd: the ground 
state of QCD matter can be analyzed in terms of the weak-coupling methods in 
QCD. Also we can count on the knowledge from condensed matter physics with 
quarks substituting for electrons. In this analogue between electrons in metal and 
quarks in quark matter, one may well anticipate that the ground state of QCD 
matter at low T should form Cooper pairs leading to colour superconductivity (CSC) 
[13, 15, 18, 19, 60, 61]. Theoretical characterization will be further elucidated in 
section 3.3 and section 6. 

There are many patterns of Cooper pairing and thus many different CSC states. 
The search for the most stable CSC state still remains unsettled except for /iq 3> Aqcd 
or TOg — > 0. It is in fact the strange quark mass ttIs that makes the problem cum- 
bersome. In the intermediate density region particularly, the Fermi surface mismatch 
(5/iq of different quark flavours is given by (5/iq ^ m'^/fj.q. When the gap energy A 
is comparable with 5/Xq, an inhomogeneous diquark condensation may have a chance 
to develop energetically. Such an inhomogeneous CSC state gives rise to a crystal 
structure with respect to A(a;), that is the crystalline CSC phase [62]. There are, at 
the same time, various candidates over the crystalline CSC phase in this intermediate 
density region and the true ground state has not been fully revealed there (see sec- 
tion 7). 

So far we have quickly looked over the key phases labelled in figure 1 and 
important points specifically picked up in figure 2. In section 3 we will proceed to 
the theoretical framework to deal with the phase transitions of quark deconfinement 
and chiral restoration, respectively. 
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3. Order parameters for the QCD phase transition 

QCD has (at least) three order parameters for quark deconfinement, chiral symmetry 
restoration and colour superconductivity as discussed below. 

3.1. Polyakov loop and quark deconfinement 

The Polyakov loop which characterizes the deconfinement transition in Euclidean 
space-time is defined as [63, 64] 



L{x) = V cxp 



-\g I dx4 A^{x^X4) 



(1) 



which is an Nc x A^c matrix in colour space. Here (3 is the inverse temperature /3 = 1/T, 
and V represents the path ordering. We will use £ to represent the traced Polyakov 
loop, 

t^—tiL. (2) 

Let us consider the centre Z(iVc) of the colour gauge group SU(A'c): elements of the 
centre commute with all SU(7Vc) elements and can be written as Zfcl with Zk = e^^^'^^-^" 
(fc = 0, 1, . . . , A^c^ 1) and 1 being an Nc x A^c unit matrix. Under a non-periodic gauge 
transformation of the following form; Vk{x) = [z^l]^''''^, the gauge fields receive a 
constant shift, 

A, ^ Al = V, [A, iig)-'d,] V^=A,-^, (3) 

so that the traced Polyakov loop transforms as i? — )• Zk(. Because A'l still keeps the 
periodicity in X/^, such a non-periodic gauge transformation still forms a symmetry of 
the gauge action. This is called centre symmetry [8, 35]. The quark action (with the 
quark field denoted by -0) explicitly breaks centre symmetry because the transformed 
field Vk{x)ip{x) does not respect the anti-periodic boundary condition any longer. 
Thus, centre symmetry is an exact symmetry only in the pure gluonic theory where 
dynamical quarks are absent or quark masses are infinitely heavy {mq — > oo). 

The expectation value of the Polyakov loop and its correlation in the pure gluonic 
theory can be written as [65, 66, 67] 

$ = {£{x)) = G-f^i'% $ == {iHx)) = e-'^^S (4) 

(^^(a;)^(y)) ^e-'3/4<'(^-2/). (5) 

Here, the constant /q (/q) independent of x is an excess free energy for a static quark 
(anti-quark) in a hot gluon medium. % Also, fqq{x — y) is an excess free energy for 
an anti-quark at x and a quark at y. § 

In the confining phase of the pure gluonic theory, the free energy of a single quark 
diverges (/q — >■ cx)) and the potential between a quark and an anti-quark increases 
linearly at long distance {fqq{r -^ oo) -^ ar with r = jcc — y|), which leads to $ ^- 
and (£^(r -^ oo)i'(O)) -^ 0. On the other hand, in the dcconfincd phase, the free 

{ Strictly speaking, the expectation value in the pure gluonic theory {(-{x)) should be defined by 
taking the limit rriq — >■ -|-oo after taking the thermodynamic limit V — >■ oo, so that it takes a real 
value. Alternatively, one may use |(^(a;))| to define the free energy of a single quark. 
§ Even when there are dynamical quarks, one may use the Polyakov loops, (4) and (5), to define the 
heavy-quark free energies. 
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Confined (Disordered) Phase Deconfined (Ordered) Piiase 



Free Energy /q = oo /q < oo 



/qq ^ err /qq -- /q + /q + a- 



Polyakov Loop (£) =0 (^) 7^ 

(r^oo) (£t(r)^(0))^0 (£t(r)^(0)) ^ |(^) p ^ Q 

Table 1. Behaviour of the expectation value and the correlation of the Polyakov 
loop in the confined and deconfined phases in the pure gluonic theory. 

energy of a single quark; is finite (/q < oo). Also the potential between a quark and 
an anti-quark is of the Yukawa type at long distance with a magnetic screening mass 
ruM [68, 69, 70], 

g-TOMr 
/qq(r -^ CX)) ^ /q + /q + « , (6) 

where a is a dimensionless constant. Note that the glueball exchange with mass of 
0{g^T) which is smaller than the electric screening scale 0{gT) dominates over the 
long-range correlation at weak coupling. Therefore, $ 7^ and {() [r -^ 00) ^(0)) 7^ 0. 

The qualitative behaviour of $ = {€) and {i)i) is summarized in table 1. We 
see that $ can nicely characterize the state of matter in such a way that i behaves 
like a magnetization in 3D classical spin systems [35, 8, 71]. The pure gluonic theory 
for Nc = 2 and 3 have been studied in lattice gauge simulations with the finite-size 
scaling analysis [72, 37]. It was shown that there is a second-order phase transition for 
Nc — 2 and a first-order phase transition for A^c = 3: For Nc = 2 the critical exponent 
is found to agree with the Z(2) Ising model in accordance with the universality. For 
Nc = 3 the Ginzburg-Landau free energy for the Polyakov loop has a cubic invariant, 
so that the first-order transition can natually be induced [8, 37, 73]. Recent studies 
of the pure gluonic theory with N^ = 4, 6, 8, 10 indicate that the transition is of 
first order for Ac > 3 and becomes stronger as A'c increases [74, 75]. This behavior is 
similar to that of the 3D Ac-state Potts model [76]. 

The idea to construct the deconfinement order parameter can be extended to a 
more general setup. Suppose that there is an arbitrary operator O written in terms 
of gauge fields that is not invariant under centre transformation. (If we take O as 
the quark propagator straight up along the imaginary-time direction, the following 
argument shall result in the standard definition of the Polyakov loop.) The expectation 
value of O in the pure gluonic theory can be decomposed by the triality projection 
[77]; (0[A])p„„ = E,'^=o'(^"[^])purc with 

(a[A])p,,c - ^ E (0[^']>Purc e-^-'^-ZA^c^ (7) 

'^ k=0 

where A'^ is defined in (3). Then, (C',i[A])pui.c with n ^ can be an order parameter 
sensitive to the spontaneous breaking of centre symmetry. A particular choice O = ipip 
[78, 79] (so-called the dual condensate) has some practical advantages [80]. 

Centre symmetry discussed so far is rigorously defined only in the pure gluonic 
system as already mentioned: In the presence of dynamical quarks, centre symmetry 
is broken explicitly, so that $ and i" always take finite values. This is analogous to 
the spin system under external magnetic field, in which the magnetization is always 
non- vanishing [81]. Nevertheless, $ may be estimated as $ = e~^'^ with a hadronic 
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mass scale (M ^ 0.8 GcV) at low T and thus the low-T phase approximately preserves 
centre symmetry. 

3.2. Chiral condensate and dynamical breaking of chiral symmetry 

In the QCD vacuum at T ~ fiB = 0, chiral symmetry is spontaneously broken, which 
is the source of hadron masses. It is a common wisdom that the chiral symmetry 
breaking is driven by the expectation value of an operator which transforms as 
(A^f, Nf) + {Nf , N[) under chiral symmetry. A simplest choice of the order parameter 
for the chiral symmetry breaking is a bilinear form called the chiral condensate, 

W) = (V^rV'l + V^lV'r), (8) 

where colour and flavour indices of the quark fields are to be summed. If the above 
chiral condensate is non-vanishing even after taking the limit of zero quark masses, 
chiral symmetry is spontaneously broken according to the pattern Q ^ % with 

g = SU(iVf)L X SU(iVf)R X U(1)b X Z(2Arf)A (9) 

n - su(iVf)v X u(i)b, (10) 

which leads to Nf — 1 massless Nambu-Goldstone bosons for Nf > 1. \\ The U(1)a 
symmetry in the classical level of the QCD Lagrangian is broken down explicitly to 
Z{2N[)a in the quantum level. Then, the U(1)a current is no longer conserved (U(1)a 
anomaly); % 

The right-hand side of the above relation is nothing but the topological charge density. 
Thus, gauge configurations with non-trivial topology are microscopically responsible 
for the U(1)a anomaly. In other words, the U(1)a current could be approximately 
conserved if the gauge configurations arc dominated by topologically trivial sectors. 
We will come back to this point later to address effective restoration of U(1)a 
symmetry in the medium [83]. 

Note that (8) is not a unique choice for the order parameter [84, 85, 86, 87]. For 
example, one may consider the following four-quark condensate, 

/v^Y(1 - 75)^^ • V5 y(1 + 75)^) = /^rA>l • ^lA'^tArV (12) 

If this is non-vanishing, the ground state breaks chiral symmetry Q to "H x Z{N{)a 
where Z(7Vf)A corresponds to a discrete axial rotation. If the bilinear condensate (8) 
is non-zero, the four-quark condensate (12) takes a finite value in general. However, 
a non-zero value of (12) does not necessarily enforce a finite value of (8). In fact, if 
Z(Aff )a symmetry is left unbroken, (8) must vanish. As long as the Dirac determinant 
in QCD is positive definite, possibility of unbroken Z(iVf )a symmetry has been ruled 
out by the exact QCD inequality [85]. However, it is not necessary the case for finite 
Hb with which the positivity of the Dirac determinant does not hold. 

II Strictly speaking, the quotient groups g' = g/[Z{N() X Z(N{)] and W = 'H/Z(Aff) need to be 
introduced to avoid double counting of the discrete centre group operation. A simplest example of 
this kind is V{N) = [S\J{N) x U(l)]/Z(Af) [82]. 

5 The anomaly relation does not necessarily guarantee the absence of the massless Nambu-Goldstone 
boson. In the Schwinger model which is an exactly solvable QCD analogue, the U(1)a problem is 
resolved by the Kogut-Susskind dipole ghosts. 
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3.3. Diquark condensate and colour superconductivity 

QCD at high baryon density shows a novel mechanism of spontaneous chiral symmetry 
breaking [88] . The fundamental degrees of freedom in the CSC phase with three colours 
and three flavours are the diquarks defined as 

where (i, j, k) are flavour indices and (a, /3, 7) are colour indices. Note that the charge 
conjugation matrix C — 17^7° is necessary to make (/?l/r be Lorentz scalar. Then, 
'y'L/R is a triplet both in colour and flavour, so that it transforms in the same way as 
the quark field V'l/r- 

Under certain gauge fixing, one may consider the expectation values of these 
operators. They are called the diquark condensates as will be discussed further in 
section 6. Instead, we can also construct an analogue of (8) in terms of the diquarks, 

{vW) + {vivK)^ (14) 

where colour and fiavour indices are summed. This is a gauge-invariant four-quark 
condensate which characterizes the spontaneous chiral symmetry breaking in CSC 
[13]. Unlike the bilinear operator in (8), the four-quark operator here keeps additional 
invariance Z(2)l x Z(2)r corresponding to the refiections; V'l ~^ ~V'l ^'^d V'r — ^ "V'rj 
which is broken down to Z(2) due to the U(1)a anomaly. 

In the CSC phase, not only chiral symmetry but also U(1)b symmetry associated 
with the baryon number conservation may be spontaneously broken. A colour singlet 
order parameter to detect such symmetry breaking can be 

ea/37eijfe((</?L/R)ai (</'L/r)/3j {'fh/K)ik)- (15) 

The six-quark operator here breaks U(1)b with its Z(6)b subgroup maintained. If 
this condensate is non-zero, there appears exactly massless Nambu-Goldstone boson 
because the baryon number symmetry is an exact symmetry in the QCD Lagrangian. 

4. Chiral phase transition at finite temperature 

The chiral phase transition at finite T with /in = has been and is being extensively 
studied by the renormalization group method near the critical point a la Ginzburg- 
Landau- Wilson and by the lattice-QCD simulations. In this section we will bricfiy 
summarize the current status of these studies. (See [39] for further details.) 

4.1. Ginzburg- Landau- Wilson analysis 

If the phase transition is of second order or of weak first order, one may write down 
the free-energy functional in terms of the order parameter field <I> as a power series 
of $/r<.. The large fiuctuation of $ near the critical point is then taken into account 
by the renormalization group method. This is called the Ginzburg-Landau- Wilson 
approach. For chiral phase transition in QCD, the relevant order parameter field is 
a iVf X A^f matrix in flavour space, ^ij ^ (V'i(l ~ 75)V'i)- Under the fiavour chiral 
rotation U(iVf)L x U(iVf)R, $ transforms as $ — > VL$Vf[. Then the Ginzburg-Landau 
free energy in three spatial dimensions {D = 3) with full U(A'f)L x U(Aff)R symmetry 
up to the quartic order in ^ij becomes [36, 47]; 

^ym- ^trV$^V$+— tr$t$+^(tr$^$)%^tr($^$)^ (16) 
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JVf ^ 2 A^f > 3 

U(1)a symmetric (cq = 0) Fluctuation-induced 1st order 1st order 
U(1)a broken (cq ^ 0) 2nd order [0(4) universality] 1st order 

Table 2. Order of the chiral phase transition conjectured from the chiral effective 
theory with niasslcss A^f flavours with and without the U(1)a anomaly. 



Effects of temperature T enter through the parameters oq, foi and 62. Note that 0.sym 
is bounded from below as long as 61 + b2/N{ > and 62 > are satisfied. The 
renormalization group analysis of (16) on the basis of the leading-order e(= 4 — D) 
expansion leads to a conclusion that there is no stable IR fixed point for Nf > v3 [36]. 
This implies that the thermal phase transition described by (16) is of the fluctuation- 
induced first order for two or more flavours. 

In QCD, however, there is U(1)a anomaly and the correct chiral symmetry is 
SU(A^f )l X SU(Aff )r, X U(1)b x Z(2A^f )a for A^f massless quarks. The lowest dimensional 
operator which breaks U(1)a symmetry explicitly while keeping the rest of chiral 
symmetry is the Kobayashi-Maskawa~'t Hooft (KMT) term [89, 90, 91, 92]; 

^anomaly = - y (dct $ + dct $^ ) . (17) 

The coefficient Cg, which is T-dependent in general, dictates the strength of U(1)a 
anomaly. In the instanton picture [91, 92], co{T = 0) is proportional to the instanton 
density ni^st, which is perturbatively evaluated as [7] 

co(T = 0)ocni„,t(p,r-0)-(^) e-«^/Sp-^ (18) 

with p being a typical instanton size. 

For A^f — 3 the KMT term becomes a cubic invariant in the order parameter. 
Hence, ri[$] = ilsym + ^anomaly Icads to the chiral phase transition of first order. For 
A'f = 2, on the other hand, the KMT term becomes a quadratic invariant. Also the 
chiral symmetry in this case is SU(2)l x SU(2)r = S0(4). Such an effective theory 
with 0(4) symmetry has a Wilson-Fisher type IR fixed point as long as the coefficient 
of the quartic term of $ is positive. Therefore, if the chiral phase transition of massless 
Nf = 2 QCD is of second order, its critical exponents would be the same as those in 
the 3D 0(4) effective theory according to the notion of universality. In Table 2 we 
summarize the Ginzburg-Landau- Wilson analysis from the chiral effective theory [36] . 

In the real world, none of quark is exactly massless: For example, TOu — (1.5- 
3.3) MeV, nid = (3.5-6.0) MeV and mg = 105^35 MeV at the renormalization scale of 
2 GeV [93] . Therefore, it is useful to draw a phase diagram by treating quark masses as 
external parameters. This is called the Columbia plot [94] as shown in figure 3 where 
the isospin degeneracy is assumed (ttJu = ttij = rriud)- The first-order chiral transition 
and the first-order deconfinement transition at finite T are indicated by the left-bottom 
region and the right-top region, respectively. The chiral and deconfinement critical 
lines, which separate the first-order and crossover regions, belong to a universality 
class of the 3D Z(2) Ising model except for special points at ?Tiud = or Wg = [95]. 

If the chiral transition is of second order for massless Nf — 2 case, the Z(2) chiral 
critical line meets the ?7iud — axis at TOs = ml"^^ (tricritical point) and changes its 
universality to 0(4) for nis > m*''' [96]. The tricritical point at rus = ml" is a Gaussian 
fixed point of the 3D (p^ model (that is, the critical dimension is not 4 but 3 at the 
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Figure 3. Schematic figure of tlie Columbia phase diagram in 3-flavour QCD at 
Hb = on the plane with the light and heavy quark masses. The U(1)a symmetry 
restoration is not taken into account. Near the left-bottom corner the chiral 
phase transition is of first order and turns to smooth crossover as mud and/or 
rris increase. The right-top corner indicates the deconfincmcnt phase transition 
in the pure gluonic dynamics. 



tricritical point), so that the critical exponents take the classical (mean- field) values 
[97] , which is confirmed in numerical studies of the chiral model [98] . 



4-2. Lattice QCD simulations 

Although the critical properties expected from the Ginzburg-Landau- Wilson analysis 
discussed above are expected to be universal, the quantities such as the critical 
temperature and the equation of state depend on the details of microscopic dynamics. 
In QCD, only a reliable method known for microscopic calculation is the lattice- 
QCD simulation in which the functional integration is carried out on the space-time 
lattice with a lattice spacing a and the lattice volume V by the method of importance 
sampling. In lattice-QCD simulations there are at least two extrapolations required 
to obtain physical results; the extrapolation to the continuum limit (a — >■ 0) and the 
extrapolation to the thermodynamic limit {V — > oo). Therefore, lattice results receive 
not only statistical errors due to the importance sampling but systematic errors due 
to the extrapolations also. 

For nearly massless fermions in QCD, there is an extra complication to reconcile 
chiral symmetry and lattice discretization; the Wilson fermion and the staggered 
fermion have been the standard ways to define light quarks on the lattice, while 
the domain-wall fermion and the overlap fermion recently proposed have more solid 
theoretical ground although the simulation costs are higher. For various applications 
of lattice-QCD simulations to the system at finite T and /ib, see a recent review [39]. 

Here we mention only two points relevant to the discussions below: (i) The 
thermal transition for physical quark masses is likely to be crossover as indicated 
by a star-symbol in figure 3. This is based on the finite-size scaling analysis using 
staggered fermion [38]. Confirmation of this result by other fermion formalisms is 
necessary, however, (ii) The (pseudo)-critical temperature Tpc with different types of 
fermions and with different lattice spacings are summarized in figure 4. In view of 
these data with error bars, we adopt a conservative estimate at present; Tpc = 150- 
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Figure 4. Determination of the pseudo-critical temperature Tpc for thermal 
QCD transition(s) from recent lattice QCD simulations. (1) 169(12)(4) MeV for 
2 + 1 flavours in the asqtad action with Nt up to 8 determined by Xm/T'^ (where 
Xm is the chiral susceptibility) [99]. (2) 192(7){4)MeV for 2 + 1 flavours in the 
p4fat3 staggered action with Nt up to 6 determined by Xm. and XL (where xl 
is the Polyakov loop susceptibility) [100]. (3) 151(3)(3)MeV and 176(3)(4)MeV 
for 2 + 1 flavours in the stout-link improved staggered action with Aft up to 10 
determined by Xm/T* and XL respectively [101]. (4) 172(7) MeV for 2 flavours 
in clover improved Wilson action with Nt up to 6 determined by xl [102]. (5) 
152(3)(3)MeV and 170(4)(3)MeV for 2 + 1 flavours in the stout-link improved 
staggered action with Nt up to 12 determined by Xm/T'^ and xl respectively 
[103]. (6) 185-195 MeV for 2 + 1 flavours in the asqtad and p4 actions with A^t 
up to 8 determined by Xm and xl [104]. (7) 174(3)(6)MeV for 2 flavours in 
the improved Wilson action with Nt up to 12 determined by Xm and xl [105]. 
(8) 171(10)(17)MeV for 2 + 1 flavours in the domain-wall action with A^t = 8 
determined by Xm./T'^ [106]. (9) 147(2)(3)MeV and 165(5)(3)MeV for 2 + 1 
flavours in the stout-link improved staggered action with Nt up to 16 determined 
by Xm/T'^ and Xa/T^ (where Xa is the strange-quark susceptibility) [107]. 



200 MeV. + It has been clarified recently that improvement of the staggered action 
with less taste-symmetry breaking favours smaller value of Tpc < 170 MeV [107, 108]. 

5. Chiral phase transition at finite baryon density 

Let us now introduce the baryon chemical potential /xb as an extra axis to the 
Columbia plot. In the so-called "standard scenario" the first-order region in the lower 
left corner of the Columbia plot is elongated with increasing ^b as written in the left 
panel of figure 5. If the physical point at ^b = is in the crossover region, there 
arises a critical chemical potential /lje so that the system shows a first-order transition 
for ^B > A*E- That is, we find a QCD critical point at (/xe,Te) on the QCD phase 
diagram in the hb-T plane. In the right of figure 5, the so-called "exotic scenario" 
is sketched, in which the size of the first-order region shrinks as /iB increases. In this 
case, if the physical point at ^b = is in the crossover region, it stays crossover for 



the physical unit, and the prescription of deflning Tpc. For crossover transition Tpc may depend on 
which susceptibility (either chiral susceptibility Xm{T) or Polyakov loop susceptibility XhiT)) are 
used, and also depend on T-dependent normalization of the susceptibilities. 
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Figure 5. Schematic evolution of the Columbia plot with increasing fiB in the 
standard scenario (left) and the exotic scenario (right). 



finite /iB, SO that there arises no critical point (at least for small /xb) in the QCD 
phase diagram in the ji-q-T plane. In general the critical surface in figure 5 can have 
more complicated structure, which may allow for several QCD critical points in the 
Hb-T plane. 

5.1. Lattice QCD at low baryon density 

In the presence of /iq ^ 0, the QCD partition function on the lattice is written as 

Z(r,M)= /[dC/]det[F(A.q)]exp-''^™(^), (19) 



where U is the matrix-valued gauge field defined on the links. The Yang-Mills action 
is denoted by (3Sym{U) with (3 = 2Nc/g'^, while the quark contribution is denoted 
by the determinant of F{fj,q) = -D(/iq) + mq with D{fj,q) being the Euclidean Dirac 
operator. The quark mass m^ is real and positive. Let ip be an eigenfunction of 
D(iiq) as D{fiq)'ipi — Xitpi with an eigenvalue A^. If /iq = then D{0) is an anti- 
Hermitean operator and thus A^ should be pure imaginary. Besides, j^ipi is an 
independent eigenfunction with an eigenvalue — A^ — A*, because of 7^-Hermiticity 
7^-0(0)7^ = ^^(0) = ^^(0). Thus, we always have complex-conjugate pairs (Ai,A*) 
and the determinant becomes real and positive; 

det[D{0) + TOq] = Y[{X, + TOq)(A* + mq) > 0. (20) 

i 

In this way, at /iq = 0, the integrand in (19) is positive definite and the importance 
sampling method works fine to evaluate the functional integral. 
For /iq 7^ the 7^-IIermiticity is replaced by 

J'D{^,qh'^DH^^^l). (21) 

Therefore dct[Z?(/iq) + 77iq] is no longer real unless /iq is pure imaginary. This requires 
us to deal with a severe cancellation of positive and negative numbers to evaluate 
the functional integral. This is the notorious sign problem whose difficulty grows 
exponentially as the lattice volume V increases. 

In the following we shall briefiy look over some of approaches in the lattice-QCD 
simulations at finite /iq. For more details, see the reviews [109, 110, 111] and references 
therein. 

• Multi-parameter reweighting method — An expectation value of an operator O at 
finite /iq can be formally rewritten in terms of an ensemble average at /iq = 0; 

(O)^, - {O ■ i?(/*q));., = 0(i?(/iq)>;, = 0, (22) 
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Figure 6. Conceptual illustration of the single-parameter reweighting method 
(Glasgow method) in the fiq direction only and the multi-parameter reweighting 
method in the both ^q and /3 (i.e. T) directions 



where i?(/iq) = detF(/iq)/dct F(0) is called the reweighting factor [112]. The 
gauge configurations generated at /iq = only occasionally sample the region 
where O ■ -R(/iq) and i?(/iq) are large (overlap problem). Also, they take complex 
values (sign problem). Therefore, such a simulation works only when /iq and V 
are small. 

To have a better overlap, the reweighting may be generalized towards not only 
the /Xq direction but also the (3 direction [113]: 

(O)^, = {O ■ i?(/.q,<5/3))^„,^^.o(i?(Mq,'5/?)>^„!^,=o. (23) 

where the multi-parameter reweighting factor is defined as R{pq, 6f3) = 
R{fiq)exp{'-S(3S) with S(3 = (3 — (3o. Here /3o should be chosen to maximize 
the overlap. The physical temperature T and f3 = IN^jg^ are implicitly related 
through the lattice spacing a. Then, the reweighting along the phase boundary 
in the /iq-T(/3) plane (as sketched in the right panel of figure 6) would have a 
better chance to probe into larger /iq. 

Taylor expansion method — The full evaluation of the reweighting factor -R(/iq) 
is not an easy task even on the computer. If /tq is small enough, the right-hand 
side of (22) can be expanded in terms of /iq/T [114, 115, 116], 



n=0 

QU 



rpn 



n 



dK 



((O.i?(/iq))^^.0(i?(Mq));'= 



(24) 



The coefficients c„ are written in terms of the quark propagator and can be 
simulated at /iq = 0. Because the expansion is based on (22), the original overlap 
problem translates into the convergence problem in higher order terms; the Taylor 
expansion makes sense for iiq/T within a radius of convergence dictated by the 
singularity closest to the origin in the complex /tq/T-plane. 



• Imaginary chemical potential method — If /iq is pure imaginary (which is here 
denoted by /Xq), (21) reduces to the 7^-Hermiticity, so that there is no sign problem 
[117, 118]. Then, it is possible to perform lattice simulations to find physical 
observables by the analytic continuation back to the real chemical potential. 






n=0 



n=0 



Mo ^ 



(25) 
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The applicability of this method is bounded by singularities or periodicity. For 
imaginary chemical potential, /iq/T plays a role of an angle variable which has 
naive periodicity by 27r. However, the partition function is a function of /iq +5A4 
and a change in /iq by 27r/7Vc can be absorbed by the centre transformation (3). 
Therefore, the actual period of /iq/T is 2tt/Nc (Roberge- Weiss (RW) periodicity) 
[119]. In the deconfinement phase at high T, especially, there is a first-order phase 
transition at a half of the RW point; /iq/T = tt/TVc ^ 1 for N^ — 3. Thus the 
method of imaginary chemical potential works up to jlq/T < 1 at best, which is 
also confirmed in model studies [120]. 

• Canonical ensemble method — In the thermodynamic limit the canonical 
ensemble with fixed particle number A'q is equivalent with the grand canonical 
ensemble with fixed chemical potential /iq (and the mean value of A^q is 
specified). To convert the grand canonical to the canonical description, the 
Fourier transform in terms of the imaginary chemical potential is necessary 
[121, 122, 123, 117, 124, 125], 



• 



2lT 



d-/- -iN^4, 



(0>A^,-y^ ^e-^''^UO)f.,=^T. (26) 

In this case, the difficulty of the sign problem is transferred to the integration 
with respect to (f) — /iq/T. This canonical approach works fine as long as the 
volume is not large, for which centre symmetry is forced to be restored for 
TVq that is a multiple of Nc [77]. It is a highly delicate procedure to take the 
correct thermodynamic limit T^ — > 00 in which not A^q but riq = Nq/V should be 
kept fixed. Especially, the order of the phase transition is sensitive to how the 
thermodynamic limit is approached [126]. 

Density of states method — There are several variants of the density of states 
method depending on the choice of a variable to rewrite the partition function. 
Here we take an example of a phase 9 of the Dirac determinant [127, 128]. (The 
plaquette P or energy E is useful as well [129, 130, 131].) The density of states 
in this case reads 

p{e) = {S{9-9{U)))^q, (27) 

where the expectation value is taken by the phase quenched simulation in which 
detT(/iq) is replaced by |detT(/iq)| to avoid the sign problem. Using p{9) one 
can express the expectation value of an operator O as 

{0),^ = ^Jd9p{9)&'{0)e, (28) 

with Z ^ Jd9 p{9) &^ and 

{0)e^^^{0-6{9-9{U))U. (29) 

Both p{9) and {0)0 are quantities calculable in the lattice simulation. The 
difficulty of the sign problem is now translated into the precise determination 
of the density of states; (21) implies that the phase quenched simulation at 
/iq 7^ with two degenerate fiavours is identical to the simulation at finite isospin 
chemical potential /xi(= Pq); 

|detF(/tq)|2 = detF{pq)dctF{-pq). (30) 
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Then, the phase quenched expectation value goes through a quaUtative change 
as /^q increases; an exotic phase with pion condensation appears for /iq > vfiT^jl 
[132]. For smah /Xq, both lattice-QCD simulation [128] and an analysis in the xPT 
[133, 134, 135] show that p{0) behaves as Gaussian, while p{Q) becomes Lorenzian 
for /iq > vfiT^jl in the xPT. In the latter case, extremely precise determination 
of piB) and (0)6/ is crucially important. So far, the density of states method is 
applied for relatively small values of /iq/T combined with the Taylor expansion 
of the fermion determinant [128, 136]. 

• Complex Langevin method — The sign problem arises from the importance 
sampling to deal with the multi-dimensional functional integral. Therefore, it may 
not appear in different quantization schemes other than the functional integral. 
The stochastic quantization [137, 138, 139, 140, 141] formulated in terms of the 
Langevin dynamics with a fictitious time is a promising candidate for such an 
alternative. If we consider a scalar field theory defined by an action S'[(/)], the 
Langevin equation reads 

5(/>(a;, s) 5S[(t)] 



ds 



'n{x,s), (31) 



where s is a fictitious time and rj a Gaussian noise; (??(x, s)) = and 
{t]{x,s)t]{x' ,s')) = 2S{x — x')S{s — s'). The expectation value is taken over the 
ri{x, s) distribution and the equilibrated value is obtained in the s — > oo limit. 
At finite /iq the action and noise are complex. The complex Langevin dynamics 
can evade the sign problem and correctly describe the system at finite density 
for some simple models, while there are also cases where it may fall into a wrong 
answer [142, 143]. Whether this method is applicable to dense QCD or not is still 
an open question. 

It would be an important milestone if the lattice-QCD simulation can show the 
existence/non-existence of the QCD critical point in p-Q^T plane. In [113, 144] the 
location of the Lee- Yang zero in the complex lattice-coupling (/3 — 2Nc/g'^) plane has 
been investigated using the multi-parameter reweighting for 2 + 1-flavour staggered 
fermions. If the thermal transition is of first order, the Lee- Yang zero nearest to the 
real axis in a finite lattice box is expected to approach to the real axis as the lattice 
volume increases, while there is no such tendency for the crossover transition. It 
was then concluded that the QCD critical point is located at {pB,T) — {pe,Te) = 
(360 ± 40MeV, 162 ± 2MeV) for physical quark masses [144]. It has been argued, 
however, in [145] that the sign problem at finite baryon density can fool the Lee- Yang 
zero near the real axis. It is also suggested to study the behaviour of a set of Lee- Yang 
zeros in the complex /3 plane to make a firm conclusion. 

If we assume that the convergence of the Taylor expansion in terms of Pq/T 
is dictated by the singularity at the critical point, the radius of convergence is an 
indicator of the location of the critical point. In [116] it was reported that the 
coefficients up to sixth order yield Te/Tpc = 0.94 ± 0.01 and pE/Tp^ = 1.8 ± 0.1 
for 2-flavour staggered fermions with mj^/mp = 0.3. In the same way, with 2 + 1- 
flavour staggered fermions, the radius of convergence is discussed in [146], though the 
results are not conclusive. The idea of the radius of convergence has been also tested 
in a chiral effective model coupled with the Polyakov loop, in which the exact location 
of the critical point and the higher-order coefficients are calculable in the mean-field 
approximation [147]. The result from the model test suggests that the relation between 
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the convergence radius of the expansion and the location of the critical point is not 
clear. 

In [148] the canonical partition function is extracted using the density of states 
method with an assumption of the Gaussian distribution of the phase 9 of the Dirac 
determinant. In this case the standard S-shape curve in the ^^~n^ plane is a signature 
of the first-order transition. The location of the critical point is estimated to be 
Te/Tpc ~ 0.76 and fiE/Tpc ~ 7.5 for 2-flavour staggered fermions with m^ ~ 770 MeV. 
In [125] the canonical ensemble method is used to study the S-shape and it was found 
Te/Tpc = 0.94 ± 0.03 and fXE/Tpc = 3.01 ± 0.12 for 3-flavour clover fermions with 
m^ ~ 700 MeV. 

As shown in figure 5 one may consider the critical surface /iB ~ I^BiiTiud, 'm-s) and 
its behaviour near /iB = to check whether the critical point is located in small /ib 
region. In particular, on the SU(3)v symmetric line (wud — "^s)? one can define the 
critical mass as a function of fiq and make a Taylor expansion, 

The sign of C2 and C4 determines whether the critical surface behaves like the standard 
scenario or the exotic scenario in figure 5. Simulations with 3-flavours of staggered 
fermions on a coarse lattice (Nt — 4) with the imaginary chemical potential method 
show C2 = —3.3(3) and C4 = —47(20) [149], which is consistent with the exotic scenario 
in figure 5. This does not, however, exclude a possibility that the critical point exists 
for large value of /iq/T. Confirmation of this result with different fermion formulations 
with smaller lattice spacing is necessary to draw a firm conclusion [150]. 

In short summary, in any method, the validity of QCD simulations at finite 
baryon density is still limited in the region /iq/T < 1 at present because of the sign 
problem. Although some of the lattice-QCD simulations suggest the existence of the 
QCD critical point in jjlb^T plane, the results are to be taken with care if it is predicted 
at large /iq/T. 

5.2. Effective U(1)a symmetry restoration 

For A'f = 3 the first-order transition at finite T is driven by the trilinear chiral 
condensates from the U(l)A-breaking KMT-type interaction with the strength cq as 
discussed in section 4.1. In the instanton calculation cq is proportional to the instanton 
density (18) which is screened by the medium at high T and /iq as 

ni„st(T,/iq) = (^)'V«-V.(p)^p-5 ,,p|_,2^2y2^2|e ^ ^^ _ ^^.^2^2 j ^ (33) 

in the one-loop calculation with the instanton size whose empirical value is p ^ 0.3 fm 
[151]. It is thus expected that the U(l)A-brcaking interaction is exponentially 
suppressed as T and/or /iq grow larger. If one performs the p-integration (which 
is not IR divergent thanks to the exponential suppression), the suppression factor 
would be T~^^ or /x" ^"^ for TVc = A^f = 3 by the dimensional reason. It is a non- 
trivial question whether the instanton-induced interaction drops by either exponential 
or power function [152]. The Columbia plot with an in-medium reduction of U(1)a 
breaking should have a smaller region for the first-order phase transition. This has 
been confirmed at /iq = [153] in the Polyakov-loop extended chiral models such as 
the PNJL model [154, 155] and the PQM model [156]. Consequently the location 
of the QCD critical point at finite /iq would be very sensitive to the strength cq 
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Figure 7. Schematic illustration of the potential shapes in the crossover, critical 
and first-order regions, respectively, from the left to the right. The i-direction at 
the critical point is tangential to the first-order phase boundary. 



[153, 157, 158]. The negative coefficients C2 and C4 in (32) might be attributed to the 
effective U(1)a restoration at larger fiq [159]; it has been found that the chiral model 
analysis with an exponential ansatz, cq oc e~ "^ ''1, leads to a scenario similar to the 
right panel of figure 5. 



5.3. QCD critical point search 

In this subsection we summarize physical consequences of the critical point under the 
assumption that it exists in the QCD phase diagram. 

Figure 7 is a schematic illustration of the shape of the effective potentials in 
the crossover and first-order transitions. If the critical point is approached along the 
first-order phase boundary, the associated potential shape looks like that of the second- 
order phase transition around a non- vanishing order parameter. The critical properties 
at the critical point can be mapped to the 3D Ising model with Z(2) symmetry with the 
(reduced) temperature t and the external field h. On the fiB-T plane in the QCD phase 
diagram the mapped t-direction is tangential to the first-order phase boundary because 
Z(2) symmetry at the critical point is locally preserved along this direction. The 
determination of the /i-direction, on the other hand, requires microscopic calculations. 

In the vicinity of the critical point along the i-direction the correlation length ^, 
the chiral condensate a and the chiral susceptibility Xa have the following scaling; 

e^t-^ a={Hj}-{{pij}o^t^, Xa^t-\ (34) 

where v ~ 0.63, /3 ~ 0.33 and 7 ~ 1.24 are known critical exponents in the 3D Ising 
model. Because chiral symmetry is explicitly broken by rriq ^ 0, the chiral condensate 
at the critical point is non- vanishing and takes a finite {•ipip)o. In the same way, along 
the /i-direction, they scale as 

C^h-''/^', a^h^", X.-h-^+^'\ (35) 

where 5 ~ 4.8. It is also possible to parametrize the singular part of the equation of 
state in a similar way as to the 3D Ising model [160]. 

The divergent correlation length S, can be interpreted physically as the vanishing 
screening mass in the (T-meson channel, which may have significant consequences in 
the HIC experiments [49, 50]. However, one should keep in mind that the pole mass 
in the a channel can never be massless since it is not directly related to ^ [161, 162]. 
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In a finite-density medium there is a mixing between the chiral condensate and 
the baryon density. Therefore the divergence in the chiral susceptibility can be seen 
also in the baryon number susceptibility, leading to enhanced fluctuation in the baryon 
number {{SN)'^). It has been argued that the general moments (cumulants) become 
[163]; 

{(SN)''), - ^K5'v)/2-3^ (36) 

where (• • -jc represents the connected piece of the correlation function. Naturally the 
higher moments have stronger divergence in the vicinity of the critical point where 
^ -^ oo. On the other hand, they require a large subtraction like {(6N)'^) — 3{{6N)^). 
It would be therefore more difficult to extract them accurately from HIC experiment 
data. 

Here we emphasize that the soft mode responsible for the critical property at the 
critical point is the density fluctuation rather than the chiral fluctuation. To see this 
point it is useful to consider the following Ginzburg-Landau expansion [164, 165]; 

F[<j,p]=-^a'-^^p' + V[a,p] (37) 

V[a, p] = aa^ + ba"^ + ca^ - ha + -/a^p + -p^ - jp. (38) 

Here p represents one of conserved charge densities such as the baryon density ne 
with an appropriate normalization. The equilibrium values of a and p arc fixed by 
dV/da = and dV/dp — 0. The dynamics is solved by the kinetic equations of motion 
dF/da = and dF/dp ~ 0, which leads to the eigen- frequencies [165], 

^- r ^'. +*T"^' ^P Ag2 ^-1+4^2^2' ^•^^' 

where Wg- and ujp are eigen-frequencies which are identified as the chiral and density 
modes, respectively. The notation Xh represents the chiral susceptibility without the 
density mixing taken into account, which diverges at the critical point. It is apparent 
from the above expressions that the soft mode at the QCD critical point is the density 
fluctuation and the a mode is a decoupled fast mode. This also explains in a natural 
way why only the screening mass of the a meson becomes vanishing at the critical 
point, while the pole mass never does as mentioned above [162]. 

5.4- First-order phase transition at high baryon density 

We now turn to the chiral phase transition at low temperature and high baryon density 
where lattice simulations are not reachable yet. In this region, only qualitative analyses 
based on various chiral effective theories are available at present. Many of these models 
predict a first-order transition in cold and dense matter; we are going to extract some 
common feature in these models based on a quasi-particle picture at /iq 7^ 0. 

Let us introduce a vacuum pressure Pvac[Afq]j which embodies the dynamical 
chiral symmetry breaking, as a function of Mq (constituent quark mass) as 

Pvac[Mq] = -a{M^ - Ml)\ (40) 

with a positive curvature parameter a. The maximum of Pvac[Afq] is realized when 
Mq = ±Mo ^ 0. (If the current quark mass ?Tiq is non-zero, a linear term in Mq, 
which favours Mq > at the maximum of Pvac[-Mq], should be present.) At zero 
temperature, as long as /Xq is smaller than Mq, nothing happens and the vacuum 
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remains empty, while rie starts arising for /Zq > Mq. The pressure from finite /iq is 
expressed as 

P^[Mq] = p dfi'^ n^ip'^) = ^ ^''^ d^^:^ (a.;2 _ M2)3/20(^:f - M^), (41) 

where nq(/j.q) represents the ferniion density in the quasi-particle approximation and 
z^=(spin)x (colour) X (flavour) is the fcrmion degeneracy. Because more particles can 
reside in the Fermi sphere for smaller mass, P^[Mq] naturally has a maximum at 
Mq = and goes to zero at Mq — /iq. 

Here -Pvac[Mq] and P^[Mq] have a peak at Mq ~ M^ and M^ = 0, respectively. 
Assuming that the total pressure is simply P[Mq] — Pvac[Mq] + P^[Mq], we can see 
that the existence of two separate peaks in P[Mq] requires [166], 

,,2 

a< -— ^4^ < -— ^ -0.076, (42) 

167r2 M§ ^ 167r2 ' ^ ' 

for 2-flavour case with u = 12. At the first-order critical point the peak at Mq = is 
just as high as the second peak at Mq ~ Mo, which determines the critical chemical 
potential, 

«-^4. (43) 

247r2 M4 "- ^ 

Once a satisfies (42), the chiral phase transition with increasing /iq at T = should 
be of first order at /Xq ~ /ic- 

The actual value of a is model-dependent: In the NJL model with two flavours 
and in the linear-cr model, a is estimated respectively as 

1 /i.A2 1 \ 

- = 



M2> ' (^^) 

—^ = 0.02 - 0.05 (linear-cr model) 

8M^ ^ ^ 

where we used the standard NJL parameters; A ~ 631 MeV, GsA^ = 2.19, and the 
resultant Mg = 336.2 MeV [10]. The uncertainty in the linear-cr model comes from the 
choice of the a meson mass. Then, in both cases, the estimated a satisfles the inequality 
(42) implying the first-order phase transition. The critical chemical potential deduced 
from (43) is, in the NJL model case, given by /ic — 1.07Mo ~ 360 MeV which is 
consistent with that obtained numerically in the NJL model. 

Let us now argue that the first-order transition obtained as above is rather 
sensitive to the choice of the model Lagrangian. Indeed, it has been known that 
the repulsive contribution to the pressure of the form +GYrii induced by a quark 
interaction of density-density type can totally wash out the first-order transition 
[167, 168, 169, 170, 171, 166]. For example, in the NJL model with the density-density 
interaction, the condition (42) is changed to 

«<16^^--3^)' (4^) 

which implies that the first-order phase transition does not arise for Gy > 0.25Gs 
[169, 170, 166]. 
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6. Formation of the diquark condensate 

Finding a ground state of quarlc matter at T « with extremely large value of /Zq 
is an interesting theoretical challenge. (We use /iq instead of ji-q throughout this 
section, for our central interest is the quark degrees of freedom) . Let us consider the 
Cooper's stability-test in quark matter [60, 61]. In the perturbative regime of QCD 
the one-gluon exchange potential is proportional to a product of the quark SU(iVc) 
charges; 

{t'')al3{t'')a'li' = Ij^ {SafjSa'fS' -Safj'Sa'fj) + —^ {SafsSa'P' + Saf3' Sa' p) -{4:6) 

The first term in the right-hand side of (46) with negative sign is anti-symmetric 
under the exchange of colour indices; a <-> a' or /3 o /3', so that a quark pair in the 
colour anti-triplet channel has attraction. The second term in the right-hand side of 
(46) with positive sign is symmetric under the same exchange of colour indices, so 
that a quark pair in the colour sextet channel has repulsion. Fermi system with two 
particles attracting each other on a sharp Fermi sphere has an instability towards the 
formation of Cooper pairs. Therefore, normal quark matter inevitably becomes colour- 
superconducting (CSC) phase with diquark condensate at asymptotic high density and 
at sufficiently low temperature [60, 61, 172, 173, 174]. 

Since quarks carry not only spin but also colour and flavour, various pairing 
patterns are possible. Hereafter we use the following notation; the colour indices a, 
(3, 7 run from 1 to 3 meaning r (red), g (green) and b (blue) in order, and in the 
same way the flavour indices i, j, k run from 1 to 3 meaning u (up), d (down) and 
s (strange) in order. We will focus our attention to quark matter below the charm 
threshold, i.e. /iq < mcharm, so that we do not take into account heavy flavours (c, b 
and t). 

Spin-zero condensate: The quark pairing with zero total spin is characterized by the 
following order parameter with 3x3 matrix structure [88]; 

{d^)ai ^ <^ap-f <^i]k{lplj C^^ ijj^k)- (47) 

The quark pair in this case is in the colour anti-symmetric (anti-triplet), flavour 
anti-symmetric (anti-triplet) and spin anti-symmetric channel to satisfy the Fermi 
statistics. The charge conjugation matrix C = 17^7" together with 7^ makes the 
above condensate a Lorentz scalar. (The effect of instantons and also quark masses 
favours the scalar condensate instead of the pseudo-scalar condensate.) If the masses 
of u, d and s quarks are all degenerate, there is an exact flavour SUv(3) symmetry. 
Then, one can always diagonalize the order parameter by the bi-unitary rotation in 
colour and flavour space to obtain {d)ia = Sia^i- We will adopt this as an ansatz even 
when u, d and s quarks are not degenerate, which is in practise a good approximation 
[175, 176]. 

Let us introduce a notation A^^s = Ai, As„ — A2, A^^j = A3 to indicate which 
quarks are involved in the Cooper pairing. If all the gaps are non-vanishing (A„d 7^ 0, 
Ads 7^ and As„ ^ 0), such a state is called the colour-flavour locked (CFL) phase [88] 
because the colour and flavour degrees of freedom are entangled with each other. For 
Ads = with other components non-vanishing, it is called the uSC phase since both 
Aud and /S.SU contain the u quark. The dSC and sSC phases are defined in the same 
way [177, 176, 178]. If only A„d is non- vanishing, it is called the 2SC (two-flavour 
superconducting) phase. In case that we refer to similar phases with other flavour 
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Table 3. Classification of the colour-superconducting phases. 

combination, we write the 2SCds, 2SCsu or 2SCud(=2SC) phase [175, 179]. When ah 
the pairing gaps are absent, the system is in the state of normal quark matter (NQM). 
We summarize these abbreviations in table 3. 

Let us now consider the symmetry breaking patterns {Q — > "H) in the CFL and 
2SC phases. * In the CFL phase for masslcss 3-flavours, one finds the following 
pattern [88, 15], 

g = SU(3)c X SU(3)l X SU(3)r xU(1)b x Z(6)a, (48) 

D U(l)em 

H- SU(3)c+l+rxZ(2)v. (49) 

'- ^ ^ 

D U(l)em+C 

Here the baryon number symmetry is broken to leave its discrete subgroup Z(2)v 
which corresponds to a reflection ip — > —ip. Also, (global) colour symmetry and chiral 
symmetry are broken simultaneously to leave their diagonal subgroup SU(3)c+l+r 
intact. Note that the electromagnetic symmetry U(l)cm associated with the electric 
charge Q = diag(2/3, —1/3, —1/3) happens to be a subgroup of the vector symmetry 
SU(3)l+r. In the CFL phase U(l)om survives as a modified symmetry, U(l)om+Cj 
which contains simultaneous electromagnetic and colour rotation. As a result, seven 
gluons and one gluon-photon mixture acquire finite Meissner mass by the Anderson- 
Higgs mechanism, while the other photon-gluon mixture remains massless. 

The 2SC phase is quite different from the CFL phase in the sense that only u 
and d quarks participate in the pairing with a chiral flavour-singlet combination. The 
symmetry breaking pattern in this case for masslcss 2-flavours with infinitely heavy 
strange quarks reads [15], 

g - SU(3)c X SU(2)l X SU(2)r x U(1)b xZ(4)a, (50) 

'- ^ ^ 

^ U(l)em 

H = SU(2)c X SU(2)l X SU(2)r x U(1)c+b ■ (51) 

3 U(l)cm+C 

Note that chiral symmetry remains intact while colour symmetry is broken, so that 
five gluons (4th to 8th) receive finite Meissner mass. In the 2SC phase, the baryon- 
number symmetry U(1)b survives as a modified symmetry, U(1)c+b7 which contains 
simultaneous electromagnetic and 8-th colour rotation. The global electromagnetism 
which is originally a combination of the baryon number symmetry and the isospin 
symmetry, survives also as a modified symmetry, U(l)cm-i-c- Therefore, the 2SC 
phase is neither a superfluid nor an electromagnetic superconductor. The modified 
symmetries discussed above remain unbroken in the uSC and dSC phases as well. 

* Similar comment on the quotient groups {Q' and W) as given to (9) and (10) is applied here. 
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The CFL phase is reminiscent of what is caUed the B phase of superfluid '^He in 
which fermionic ^He atoms have a spin-triplet (S* = 1) and p-wave (L = 1) pairing. The 
B phase is a state in which the orbital angular momentum L and the spin S are locked 
so that the symmetry breaking pattern is S0(3)s x S0(3)l x U(1)0 — > S0(3)s+l, 
where \]{1)^ is a symmetry corresponding to U(1)b in quark matter. 

Spin- one condensate: There is also a possibility of flavour symmetric Cooper pairs 
such as the pairing within the same flavour [f72]. In this case s-wave Cooper pair 
must be spin triplet (J = L-|-S' = 0-|-l = l)to satisfy the Fermi statistics. 

Various forms of the spin-one CSC states have been studied so far [179, 180, 181, 
182]. Since the Cooper pair is triplet in both colour and spin, the order parameter 
becomes a complex 3x3 matrix, which is similar to the situation in the CFL phase 
where the Cooper pair is triplet in both colour and flavour. The explicit form would 
be 

{d^)c^ ^ e„^^(4C7V7): (52) 

where i refers to not the flavour but the Lorentz index. The colour-spin locked (CSL) 
ansatz, {d)ia c>c (5iQ, leads to a symmetry breaking pattern which is similar to the CFL 
phase in CSC and to the B phase in superfluid "^He; 

SU(3)c X SU(2)j X U(l)e„, ^ SU(2)c+j. (53) 

On the other hand, if we take an ansatz, {d)ai oc SasiSn + 'i5i2), we have a symmetry 
breaking pattern, 

SU(3)c X SU(2)j X U(l)em ^ SU(2)c X Uj+em X U(l)em+C, (54) 

which is analogous to the A phase of superfluid "^He in which the symmetry breaking 
patter is; S0(3)s x S0(3)l x U(1)0 -^ U(l)s^ x U(1)l^+0. Other pairing patterns 
such as the polar phase, the planar phase and so on have been also investigated [182]. 

Weak- coupling results: In the asymptotically high-density region, the running 
coupling constant of the strong interaction becomes small enough to justify the weak- 
coupling calculations, so that the diquark condensate as well as the gap energy of 
quarks can be estimated from the first-principle QCD calculation. The driving force 
of the diquark condensate is the interaction between quarks due to gluon exchange. 
The chromoelectric part of the interaction in high-density quark matter is screened 
by the Debye mass m^ = ^/NfJ{2^^^ g ^q, while the choromomagnetic part of the 
interaction is screened only dynamically by Landau damping [183]. Therefore, instead 
of having the standard BCS form, A ^ ^q cxp(— const. /g^), the gap energy at the 
Fermi surface is parametrically enhanced due to large forward scattering between 
quarks [183, 184, 185, 186, 180]; 

Ap = 512.4 (^Ay'^V('^^+^)/« ^^a.^a.yi/2 g-,^^ exp(^-^), (55) 

where Ai = ai = 1 and A2 = 02 = for the 2SC phase, while Ai = 4 and A2 = 1 
with ai = 1/3 and 02 — 2/3 in the CFL phase. In the exponent of (55), unusual 
dependence on the coupling g~^ instead of the standard dependence g~^ enhances the 
gap substantially in the weak coupling regime both for 2SC and CFL phases. We see 
that A^sc jjj ^Yie 2SC phase is greater than that in the CFL phase as A^^C = 2^/3 A^^^, 
which holds not only in the weak-coupling QCD calculation but also in effective model 
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Figure 8. Fermi surface mismatch in unpaired normal quark matter. 



calculations. Similarly to the BCS theory, the gap energy at zero temperature and the 
critical temperature are related to each other in the mean-field approximation [180], 

rc = (AiU^^)i/2 ^Ap. (56) 

TT 

Note that the factor (A°^A2^) in the right-hand side is cancelled by the same factor in 
(55), so that A^'^c > ^cfl ^^^^ ^^^ necessary imply T^^^ > T^^^. 

6.1. Neutrality in electric and colour charges 

For the bulk system to be stable, the total electric and colour charges must be zero 
due to the Gauss law [187, 188, 189]. In some effective models which do not have 
gauge fields as explicit degrees of freedom, the charge neutrality must be imposed by 
electric and colour chemical potentials, ^ai — f'q + lJ-Q{Q)ii + ^3(73)00 -I- Hs{^Tg,)aa- 
Their explicit forms read 



fJ^ru 


= Mq - Ia^o + ^1^3 + |a*8, 


tLgd 


= Mq + |Mo - 5M3 + |M8, 


l^bs 


= Mq+ |Mc- ImS, 


fJ-rd 


= Mq+ |Mc + 5M3 + 5M8, 


l^gu 


2 111 
— /^q — 3/Uc — 2^3 ' 3^8, 


fJ-rs - 


= A*q+ |Mc+ |Ai3+ |M8, 


t^bu 


= Mq - sMc - 3M8, 


fJ-gs 


= Aiq+ |Mc- 5^3 + |M8, 


t^bd 


= Mq + 5A*o - Ia^s, 







(57) 



where the electron chemical potential jj,e satisfies /ic — —^Q due to /3-equilibrium. 
These constraints add some varieties in the physics of CSC because one needs to 
consider the Cooper pairing between quarks with mismatched Fermi surfaces. 

For unpaired quark matter with ?Tiud = and mg 7^ 0, the number of strange 
quarks is less than that of other quarks. There arise electrons in the system to 
maintain charge neutrality accordingly. The neutrality conditions obtained by the 
variation of the free energy, ^l^unpairod/^Aic = 9f^unpairod/<9/i3 == dfl^npaircd/dfis = 0, 
lead to //3 = /ig = and fi^ ~ '7ig/(4/iq). With the free dispersion relation at large 
chemical potential, eud(p) = |p| and es{p) = \/\p\'^ + ml ~ |p| -t- mg/(2/iq), the Fermi 
momenta have mismatches as 

„ ""^l d ™? , 5m^ ,^„, 

^^='^<'-6Mq' ^^ = ^^+12;^' '- = ^--W.- ('') 

This is illustrated in figure 8. 

At high density where the expansions in terms of ms//iq and A/^q are valid, 
model-independent conclusions can be drawn for the CSC quark matter at T = 
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[187, 188]. In the CFL phase the modified electric charge Q associated with the 
symmetry U(l)om+c is conserved. Also, all the Cooper pairs have Q = 0, so that the 
CFL quark matter at T = is a Q-insulator [190]. Thus, the CFL free energy ficFL is 
independent of the corresponding chemical potential except for a small contribution 
proportional to /i^ from Q-violating electrons. Under the neutrality conditions, 
dricFh/d^c = dO.cFh/dfJ'S = di^CFh/dfis = 0, we obtain 

1 Tfl 

M3 = Mc, A^s = 2^^°^ 2"^' Mo = 0. (59) 

It is an important property that the CFL phase (A„d — ^ds — ^su with /io = 0) 
is rigid against the Fermi surface mismatch [190] as long as the CFL phase is in the 
stable region; A > TOg/(2/Xq). We will discuss later what happens once this stability 
condition is violated in section 7.2. 

In the 2SC phase the colour and charge neutralities lead to 

2 
Ul 

^J■3 = ^J■8 = 0, Mo = TT^, (60) 

so that electrons are required to be present with even larger chemical potential than 
the normal quark matter. Also, the 2SC phase contains ungappcd Q-carrying modes, 
and is therefore a Q-conductor. 

6.2. Ginzhurg- Landau approach 

To identify the phase structure of CSC near second-order or weak first-order transitions 
at finite T, the Ginzburg-Landau approach is quite useful [61, 191]. If we assume that 
the gap energy is small compared to the critical temperature, the Ginzburg-Landau 
free energy with a power series expansion in terms of {d)ia up to the quartic order 
reads [177], 

Vl = atT:dd) + Pi{tidSf + (32 tridd'^f. 



■■^{\dcu\ 



|2 , IJ |2\ , V V^/i . i2 



do.d?) + lY.i\do.d? + \do.s? - 2|d„„|2), (61) 



where tr is taken in flavour space. For Trig ^ pL^ — Q the transition from the CFL 
phase to normal quark matter is driven by the parameter a changing the sign from 
negative to positive. The fourth term with e takes care of the asymmetry between (u, 
d) and s introduced by mg ^ 0, while the fifth term with 77 represents the asymmetry 
between u and (d, s) from the charge neutrality effect by ji^ 7^ 0. The effect of colour 
neutrality is negligible near T^. Under the diagonal ansatz on the gap matrix, the free 
energy simplifies into 

n^a' (AL + A^^ + AL) - e AL - r; A^^ 

+ /3i ( AL + A^^ + AL) ' + /32 ( AL + A^,, + A^) (62) 

with a' — a + e -\- rj/'i. In the weak-coupling analysis one can show that Pi = fi2 ^ P 
and a = aQ{T — Tc)/Tc, where ap = 2/iq/7r^ and Tc is the critical temperature defined 
for e — rj — 0. In the weak-coupling limit of QCD with mg <S^ /Xq, it is shown that 
e ~ 2r; ~ 2q:o(5 with a parameter S ~ (rn^/SfiV) ln(/iq/Tc). Then the gap energies near 
Tc behave as 

AL(7^)^g('^ + ^^), AL(r) = AL-^, ALm = AL-^.(63) 
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Figure 9. Ordering of the averaged Fermi surfaces for respective Cooper pairs. 

The inequalities, Aud{T) > Ads{T) > As„(T), imply sequential melting of CSC as T 
increases at high baryon density, i.e. CFL — > dSC — > 2SC -> NQM. This scenario of 
melting pattern has been confirmed by explicit calculations in the NJL model [178]. 

The ordering A„d > Ads > Asu in (63) can be understood in an intuitive manner. 
Because T is close to Tc , the Fermi surface is blurred and the mismatch of the Fermi 
momentum no longer imposes a pressure to take the pairing apart. The magnitude of 
the gap energy is then determined by the density of states. For example, the magnitude 
of the density of states for the pairing Aud is related to the size of the "averaged" 
Fermi surface over u and d quarks. The Fermi surface ordering as shown in figure 8 
gives rise to the ordering of the averaged Fermi surfaces as shown in figure 9, which 
explains the ordering A^^d > Ads > Asu near T^. 

In some model calculations in the literature, uSC instead of dSC appears even for 
extremely large /iq. This is caused by an artifact of the UV cutoff A in the model: If 
A is not chosen to be sufficiently large, e and rj suffer from significant cutoff-artifact, 
so that rj may even have an opposite sign from the QCD prediction [178]. On the 
other hand, as fj,q becomes small and approaches to mg, the expansion in terms of 
rris/ fj.q adopted for the computation of e and i] in (61) is no longer valid, and the other 
melting pattern, CFL — ^ uSC -^ 2SC — )■ NQM, indeed becomes possible [178]. 

6.3. Quark-hadron continuity and U(1)a anomaly 

The Ginzburg-Landau approach can be extended to incorporate the competition 
between the diquark condensate and the chiral condensate [192]. In this case it is 
necessary to distinguish the right-handed and left-handed diquark condensates dj^ and 
d^ as given in (13) to construct properly the Ginzburg-Landau free energy with the 
symmetry Q in (48). The basic transformation properties of the fields $, d^ and dp, 
are 

^^V^^VI d^^V^d^V^, d^^V^d^V^^ (64) 

where Vl, ^r and Vc correspond to U(3)l, U(3)r and SU(3)c rotations, respectively. 
In a simple case where all u, d and s quarks are massless, the chiral part has a 
standard expansion in the same form as (16); 

f}$ = ^tr$t$ + ^(tr$t$)2^^tr($t$)2_|l(det$ + dct$t). (65) 

The term with the coefficient cq originates from the axial anomaly. 
The diquark free energy up to the quartic order reads 

rid = atT{d^dl + d^dl^) + Pi {tr d^dl)^ + (tr d^d\^y 
+ /3jtr(dL4)'+tr(dX)' 
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Figure 10. Schematic figure on the reahzation of the quark-hadron continuity 
induced by the U(1)a anomaly. (Left) First-order phase boundary near T = 
without the anomaly-induced ctA'^ term. (Right) Smooth crossover and the 
appearance of a critical point due to a large ctA^ term with 7 > in (68). 

-f- /33 tr(dj^4'^L^R) + Pi *!■ '^l4 ti" ^rC^r- (66) 

The transition from the CFL phase to normal quark matter is driven by the parameter 
a changing the sign from negative to positive. UnHke det $ in (65), terms such as 
detc^L and detdj^ are not allowed in (66) since they break U(1)b. 

The coupling between the diquark and the chiral condensates has the following 
general form up to the quartic order; 

+ Ai tr(dL4**^ + -^R^*^*) + ^2 tr(dL4 + d^d\^) tr$^$ 

H- A3 (det $ tr ^l^*"^ + det $^ tr ^^4*') ■ (67) 

The term with the coefficient 71 originates from the axial anomaly. 

In the massless 3-flavour limit, one may assume $ = diag(cr, a, a) and d^ = —d^ = 
diag(A, A, A) and then the sum of above three pieces amounts to 

.3 , ^4\ , /^"a2 , ^a4\ .a2_ , xa2_2 



"^^ ^ U^ 3^ + r ) + l^f ^ + 5 A^ j - 7 AV + A AV^ (68) 

Here the a^ and the A^cr terms originate from the axial anomaly. From the Fierz 
transform of the anomaly-induced KMT interaction in the quark level, the coefficients 
c and 7 turn out to have the same sign [192]. Furthermore, the sign of c and 
7 can be taken to be positive without loss of generality because of the relation, 
03F(cr, A; c, 7) = J73f(— (T, A; — c, — 7), which is a situation analougous to the quark 
mass term. 

It should be noted that the A^cr term with positive coefficient 7 favours the 
coexistence of the chiral and diquark condensates. Also, this term is linear in a and 
behaves as if it is an explicit chiral symmetry breaking term. In contrast, A^cr^ 
term with positive coefficient A (the positive sign is supported by the weak-coupling 
calculation and in the NJL model) disfavours the coexistence. Therefore, if the effect 
of 7 is sufficiently strong and A is sufficiently large, the first-order phase boundary, 
which normally separates the chiral symmetry breaking phase {a ^ 0) and the CFL 
phase (cr = and A ^ 0), can be smeared out. In such a case the two phases are 
connected smoothly to each other and a critical point associated with this crossover 
appears as shown in figure 10. 

The smooth crossover of the CFL phase and the Nambu-Goldstone (hadronic) 
phase may have a close connection to the idea of quark-hadron continuity [175, 193]. 
Since the axial anomaly tends to enhance (reduce) the first-order transition through 
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the term proportional to c (7) in (68), it is a dynamical issue whether the phase 
diagram as sketched in figure 10 is realized or not in the real world [194, 195, 196]. 

6.4- Collective excitations 

An interesting and related question associated with the U(1)a anomaly is the fate of 
collective excitations. Using the chiral effective Lagrangian approach [197], one finds 
the dispersion relations for the pions and kaons in the CFL phase with Trig 7^ and 
/Uc 7^0 [198, 199], 



:(P) :-±Aie + y^«V+M2±, 



eK± (p) = ±Aic T ^ + V «'P' + ^l± ' 

^Ko ip) ^-^ + Jv^p^ + Mf,,, (69) 

where w^ = 1/3 at high density. The CFL-meson masses are given by 

M^± = a(TOu + md)ms + x(™u + md), 

M'^± == airriu + ms)md + x(™u + m^), 

M^-o = a{md + ms)mu + ximd + m-s) . (70) 

Here a = 3AV(7r^/2) with /^ = (21 - 81n2)/i2/(367r2) at high density and x 
parametrizes the contribution of the U(1)a anomaly which generates (iptp) a-nd 
therefore contributes to the CFL-meson masses. 

In the absence of the U(l)A-breaking term (x = 0), the energies for K^ and K'^ 
become negative and kaon condensation occurs for nig > m^'^A'^'^ with m being either 
TUu or TTid. In particular, the electron contribution to the thermodynamic potential in 
the CFL phase favours the K" condensation (the CFL-K^ phase) [200, 201, 202, 203]. 
The phase structure with the CFL-K'^ state and its variants have been also investigated 
in the NJL-type model [204, 205]. The onset of the K'^ condensation depends on the 
strength of the U(1)a anomaly x as is evident from (70). 

In view of (70) the meson masses have the ordering M^± > Mk± ~ M^o for 
TOg ^ niu ~ rrid and x ~ 0, which is inverse of the ordinary ordering in the vacuum 
[198]. This is, however, natural from the diquark picture as already implied by the 
order parameter (14) in which CFL-cr meson consists of two diquarks, qqqq. The 
Nambu-Goldstone bosons are accordingly composed from qqqq; CFL-tt"*" contains a ds 
diquark that transforms like u quark and an su diquark like d quark, while CFL-_ftr+ 
a ds diquark and a ud diquark like s quark. Therefore CFL-ii'+ has a d quark instead 
of an s quark as compared to CFL-7r+ and thus it becomes lighter than CFL-7r+. The 
effect of X 7^ in (70), on the other hand, favours the standard ordering. This is 
because x arises from the instanton interaction and induces a mixture of qqqq and 
qq, which is embodied in the A^cr term in (68). If the quark-hadron continuity is 
realized, not only the pseudo-scalar mesons but also the vector mesons and fermions 
would obey the spectral continuity. For example, the continuity of flavour-octet vector 
mesons at low density and the colour-octet gluons at high density together with the 
fate of the flavour-singlet vector meson have been investigated in the in-medium QCD 
sum rules [206] . 

In order to clarify the existence of another critical point associated with the 
quark-hadron continuity and the formation of if" condensation in the CFL phase, it 
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Figure 11. Particle-hole pair near the Fermi surface generating an 

inhomogeneous chiral condensate with a net momentum, the magnitude of which 
is given by 2fj,q. 



is demanded as a theory task to quantify how much the U(1)a symmetry is effectively 
restored at finite T and/or fj,q. 

7. Inhomogeneous states 

In the intermediate density regions of the QCD phase diagram, the ground state may 
have inhomogeneity with respect to the condensates. The conventional 7r~ and tt^ 
condensations induced by the p-wave pion-nucleon interaction in nuclear matter and 
in neutron matter are well-known examples [207]. In this section, we will address 
some of the proposed inhomogeneous phases associated with chiral transition and 
with colour superconductivity. 

7.1. Chiral- density waves 

In various materials typically at low dimensions, the charge-density wave (CDW) 
[208] and the spin-density wave (SDW) [209] are realized as ground states. Their 
counterpart in quark matter is the chiral-density waves: Such a possibility has been 
discussed in the large N^ limit of QCD [210, 211] with a spatially modulated chiral 
condensate, 

{{p{x)ip{x)) = acos{2q ■ x), (71) 

with |q| — fiq. Figure 11 is a schematic picture of the pairing of a quark and a 
quark-hole which makes an inhomogeneous condensate. The magnitude of the net 
momentum is given by 2/iq. This pairing is favoured by the forward scattering with 
singular interaction ^ l/q^ at high density. The conclusion of [211] is that this chiral- 
density wave state can be realized only for Nc > lOOOiVf . The reason why large A^c is 
required is that the gluon interaction ~ 1/q^ is IR screened (Thomas-Fermi screening 
and Landau damping) by the quark loops of 0{1/Nc). 

In [212] a different ansatz in 2-flavour case has been investigated; 

{^p{x)ij(x)) - i{i;{x)i-f^T3tP{x)) ^ae'^'"''. (72) 

Here T3 is a third component of the flavour Pauli matrices. This type of condensate is 
known to be realized as the chiral spiral in (1 + l)-dinicnsional chiral models at finite 
density [213]. Recently the chiral spiral in QCD is analyzed with confinement effects 
taken into account [214]. The confining interaction ~ l/{q^)'^ is more IR singular 
than the perturbative interaction, so that the pairing like drawn in figure 11 may be 
realized for smaller N^. than discussed in [211]. 
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In addition to the IR singular forward scattering, there is another source to 
produce a spatial modulation of the chiral condensate. To see this, the Ginzburg- 
Landau approach is again useful. The general terms of the free energy up to the sixth 
order including derivatives are [215] 

VI = C2M^ + CiM'^ + 4(VM)2 _^ ^^^6 _^ c'^iyMfM"^ + Ce(AM)2, (73) 

where M is an order parameter for the chiral phase transition (either the chiral 
condensate or the constituent quark mass in the chiral limit). Note that the magnitude 
of M and its spatial variation are assumed to be the same order. If C4 is negative, 
inhomogeneous M tends to be favoured. However, the sign of c\ cannot be determined 
by symmetry argument alone. In the mean-field treatment of the 2-flavour NJL model, 
the expansion turns out to be 



"2 ,,,2 , «4r,,4 , m./n21 , "6 



1 



n = -^M^ + -i M-* + (VM)' + -^ M" + 5(VM)'Ar + -{/^MY , (74) 
2 4 6 L 2 J 

where a2, 0:4 and ag Eire all calculable within the model. In particular, the sign of 

Ui controls whether the inhomogeneous state occurs: If it is negative, the first-order 

transition is driven and simultaneously a spatial modulation develops. Therefore the 

first-order phase boundary is naturally surrounded by the inhomogeneous state in the 

QCD phase diagram as illustrated in figure 1. Essentially the same phenomena takes 

place in the chiral Gross-Neveu model in (1 + l)-dimensions [213]. 

Importance of inhomogeneous phases in the QCD phase diagram has been 

revisited recently and physics implications have not been fully unveiled. Especially it 

is an urgent but unanswered question how the inhomogeneous condensate could affect 

the region around, if any, the QCD critical point. 

7.2. Implications from chroniomagnetic instability 

A similar inhomogeneity arises not only in the chiral condensate but also in the diquark 
condensate A (a;). The CFL phase is rigid and /S.ud — ^ds — ^su — A as long as 
(5/iq = rn^/ ^q < 2 A. Once the Fermi surface mismatch J/iq exceeds this bound, a new 
form of colour superconductivity appears, which is a counterpart of what is known 
as the Sarma state in condensed matter physics. The Sarma state is, however, not 
stable as it is, and it has been argued that a number constraint such as the charge 
neutrality condition may help the stabilization [216]. In the QCD context [217] such 
superconducting states are called the gapless 2SC (g2SC) phase [218, 219] and gapless 
CFL (gCFL) phase [220, 221] in the 2-flavour and 3-flavour cases, respectively. 

It has been found, however, that the gapless superconducting states suffer from 
another instability problem. In [222, 223, 224, 225, 226] the Meissner (magnetic 
screening) masses in the g2SC and gCFL phases have been calculated and turned out 
to be imaginary. That is, there appear negative eigenvalues from the mass-squared 
matrix in colour space. 



3 £f dAfOAf 



(75) 

for transverse gluons Af. This fact is commonly referred to as the chromomagnetic 
instability. The physics implication of the chromomagnetic instability can be nicely 
articulated by using the Ginzburg-Landau description. A gauged kinetic term can be 
added to the free energy in a form of 

0[A, A] - r!o[A] - K''\d, - i5A0A]t«[(a' - igA')^]\ (76) 
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Figure 12. An example of the mean-field model phase diagram with regions 
suffering from chromomagnetic instability [227]. At high density the ordering is 
the CFL, dSC and 2SC phases from the bottom to the top as explained in figure 9. 
At lower density the uSC phase is favoured because the Fermi surface mismatch 
between d and s quarks is the largest and A^j^ tends to melt first when the critical 
T is not large. The diagram is drawn with sufficiently large UV cutoff A = 1 GeV 
and A ~ 40 MeV is chosen at /iq = 500 MeV, so that the results are robust against 
cutoff artifacts; otherwise unphysical structures would easily enter. 



from which one can derive the Meissner mass-squared as 



(m^^"'3 



'M 



2K'^^(T")ca(T'3)fcrfA^A'^. We have k"^ < leading to {mlj)^^ < in the g2SC phase 



,22 



< leading to {m\,j)'^'^ — ('ti^)^^ < whose onset is shghtly de- 



and also k^^ = k" 

layed after the gapless onset. The remaining three gluons are unscreened. On the 
other hand, in the gCFL phase, all eight eigenvalues of the mass-squared matrix can 
potentially become negative [224, 226, 227]. The unstable regions are mapped out in 
the NJL model [227] onto the phase diagram, which is presented in figure 12. The 
shaded regions are unstable with respect to gluons with their colour labelled aside. 

Now let us consider a modulation of the diquark condensate in the simplest form 
of the plane-wave type. This is written as 

A(x) = |A|e'^°'J°-", (77) 

which is known as the (coloured) Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state 
[228, 229]. Then, the stability with respect to growth of q" 7^ is deduced from 
the potential curvature; 



1 



3 



on 



3 S dqrdqf 



(mlf) 



,Q/3 



(78) 



Therefore the chromomagnetic instability leads to the FFLO state. In the g2SC 
case with k^^ < a non-coloured {a = 0) FFLO state may be enough to cure the 
chromomagnetic instability [230] (see also [231] for results in opposition), but it is a 
subtle problem if the non-coloured FFLO can stabilize in all colour channels of gluons 
for 3- flavour quark matter [232]. 

Generally the instability analysis tells us a tendency towards the destination at 
best, but cannot go into the ground state. There are several candidates proposed 
so far, among which the free energy should be compared to sort out the most stable 
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Figure 13. Pairing with a Fermi surface mismatch by Sfiq. Left) BCS pairing 
with vanishing net momentum. One particle is lifted up from the inner Fermi 
surface by S/iq, which costs an energy. Right) FFLO pairing between particles 
having momenta p + q and — p + q with a net momentum 2q. The energy costs 
are reduced around the dotted rings depicted in the figure. 

ground state, which is still an open question. Below some of proposals are briefly 
reviewed. 

Crystalline colour superconductivity: It is intuitively understandable that the plane- 
wave FFLO state is a likely alternative of the ground state if the Fermi surface 
mismatch grows large. Naturally one can anticipate that \q\ should be of order of 
6fj,q and actually |q| = 1.2(5/Zq is found in the weak-coupling analysis [62]. The energy 
gain comes from the ring regions on the Fermi surface as indicated in figure 13. Thus, 
it is likely that a superposition of multiple q's could further decrease the energy by 
covering the Fermi surface with rings. The general ansatz is then, 

{u{x)d{x)) ^ A^^ Y, e^'-Ja- (79) 

for two flavours and, in addition to this, two more condensates, 

(s(x)u(x)) ^ A,„ Y. e"'"^"' {d{x)s{x)) ^ Ad, Y e""""' (80) 

q?e{'32} 9ie{qi} 

for three flavours [233]. Such a crystal structure has been investigated by means of 
the Ginzburg-Landau approach [234] with an approximation that A^^ ~ Asu and 
{d{x)s{x)) ~ because the Fermi surface mismatch between d, u and s quarks is all 
the same and the d-s pairing is least favoured (see figure 9). Then, qg and q"^ are 
optimized so that the free energy becomes smaller. In the 2-fiavour case the ground 
state takes the face-centred-cubic (FCC) structure with eight ^3 vectors. In the 3- 
flavour case, a most stable structure is the 2Cube45z that forms two cubic sets by 
^2 and ^3 with a relative rotation by 45° around the z axis. Another most stable 
structure is the CubeX where two rectangles by Qj and qg cross in the X shape to 
form a cube. It is desirable to study the crystallography in microscopic models: A 
recent step towards such direction can be seen e.g. in [235]. 

Gluonic phase: The most straightforward interpretation of the chromomagnetic 
instability would be the gluonic phase [236] in which gauge fields have a finite 
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expectation value. If only one component of gluons condenses, it is equivalent with 
the single plane- wave FFLO state. In fact, an ansatz, 

^^ = ^g{Al), M3=.g(^g>, 9{Al)^0, (81) 

is considered and named the gluonic cylindrical phase in [237]. Because only one 
spatial gluon condensate is involved, this state can be mapped to a FFLO-typc pairing 
with a coloured phase factor. On the other hand, multi-gluon condensation cannot 
be transformed into a single plane-wave FFLO state. In the 2-flavour model the free 
energy of the preferred gluonic phase has turned out to be smaller than that of the 
FFLO state in a wide parameter region [238] . Such an example is the gluonic colour- 
spin locked (GCSL) phase that is defined by the following ansatz, 

M8 = ^.9(A«), g{Al)^g{A^:)^Q, (82) 

which is free from the chromomagnetic instability [239]. It is then found that these 
gluonic phases have a smaller free energy than the single plane- wave FFLO state and 
the unstable g2SC phase. In most regions the GCSL is the most stable apart from the 
vicinity of the first-order phase transition to normal quark matter where the gluonic 
cylindrical phase is energetically favoured. It has not been understood, however, how 
to generalize the description of the gluonic phase to the 3-flavour case. Moreover, an 
extension of the gluonic condensation with spatial inhomogcneity may further decrease 
the free energy [240]. 

Meson supercurrent state: As we have already emphasized, the CFL phase 
spontaneously breaks chiral symmetry, where low-energy excitations are described 
by a chiral effective Lagrangian expressed in terms of colour-singlet modes. Then 
it is possible to formulate the instability by using the chiral effective Lagrangian 
[241]. The origin of the instability should be common, but it is no longer the 
"chromomagnetic" instability since gluons do not appear explicitly but all the physical 
degrees of freedom are Nambu-Goldstone bosons. Then, instead of gluons, one may 
expect the condensation of vector fields given by mesons, that is, the meson currents. 
Such a destination of the ground state is called the meson supercurrent state. The 
description looks different at a glance from the chromomagnetic instability and the 
single plane-wave FFLO state, but the underlying physics must be closely related. 
An advantage of using the chiral effective Lagrangian is that the inclusion of the K'^ 
condensate is straightforward, which is complicated in microscopic models. Such an 
interpretation as the supercurrent generation is also proposed in [242] not relying on 
the chiral effective Lagrangian but in terms of phase fluctuations around the diquark 
condensate. 

Mixed phase and phase separation: All the states as we have seen so far should 
compete with a rather conventional possibility, i.e. the mixed phase. Neutrality 
conditions with respect to gauge charge enforce the Fermi surface mismatch, but these 
conditions may be relaxed by a formation of a mixture of CSC and NQM regions as 
discussed in the 2-flavour case in [243]. 

To clarify the mixed phase structure, however, the precise determination of the 
surface tension is indispensable. The balance between the surface tension and the 
Coulomb energy fixes the typical domain size of the mixed phase structure. Naturally, 
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Figure 14. Schematic phase diagram of large- A^c QCD. The pressure is 0{N^) in 
the hadronic phase, 0{Nc) in the quarkyonic matter, and 0{N^) in the deconfined 
phase. 



the larger the surface tension is, the larger the favoured domain size grows, and 
eventually the mixed phase should be rather regarded as the phase separation [244]. 
It is not easy to extract the information on reliable value of the surface tension in the 
intermediate density region, however. So far, the possibility of the mixed phase in 
3-flavour CSC phase has not been studied seriously. 

8. Suggestions from QCD-like theories 

In the research towards the phase diagram of dense QCD, some knowledge from QCD- 
like theories would provide us with a useful hint to attack the QCD problem. There 
are many such attempts and it is impossible to cover all of them in this article. Here, 
we discuss some selected topics which are highly relevant to the understanding of the 
dense-QCD phase diagram. 

8.1. Quarkyonic matter at large Nc 

The novel QCD phase structure in the large A^c limit has been recently proposed [58] 
as schematically shown in figure 14. When /iq is smaller than the threshold of the 
constituent quark mass Afq ~ M-q/N^ ^ 0{N^), we have the hadronic phase with zero 
baryon density at low temperature. As the temperature is increased, there appears 
the first-order deconfinement transition at Tj ~ Aqcd at which the number of degrees 
of freedom and the pressure jump discontinuously from 0{N^) to 0{N^). Since quark 
loops are suppressed by l/Nc as compared to gluon contributions, Tj is independent 
of fj,q in this region as shown in figure 14. 

When ^q becomes greater than Afq, a non-zero baryon density is turned on. The 
pressure associated with this threshold at /iq ~ Mq changes discontinuously from 
0{N^) to 0{Nc). Also Td does not change with increasing /iq as long as /iq ^ 0(iV°), 
so that quarks are still confined in the right-bottom region (/iq > Afq and T < Td) 
in figure 14. The confining phase with the pressure of 0{Nc) is called the quarkyonic 
matter [58]. 

If we describe the quarkyonic matter as a weakly interacting quark system, the 
pressure of 0{Nc) is a natural consequence, but it is difficult to reconcile it with 
the confining feature. If we describe the quarkyonic matter as a baryonic system, 
it must be a strongly interacting matter where the pressure is dominated by baryon 
interactions of 0{Nc) rather than the kinetic pressure of 0(1/A'c). A possible idea 
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Figure 15. Intuitive picture of the Femii sphere which accommodates quarkyonic 
matter — quark Fermi sea and baryonic Fermi surface. 



to unify these two descriptions proposed in [58] is illustrated in figure 15; the quarks 
deep inside the Fermi sphere are weakly interacting, because it is hard to excite these 
quarks above the Fermi sea due to Pauli blocking. On the other hand, the quarks near 
the Fermi surface with a shell-width ~ Aqcd are not affected so much from the Pauli 
blocking and can interact strongly through IR singular gluons at large N^. Thus, the 
bulk thermodynamics such as the pressure, entropy and so on are dominated by the 
quarks inside the Fermi sphere, while the physical excitations on top of the Fermi 
surface are dominated by colour-singlet mesons and baryons. 

Whether the remnant of quarkyonic matter at large Nc remains in the QCD 
phase diagram at A'c = 3 is an open question. Also, how the chiral transition takes 
place inside quarkyonic matter is an important problem to be studied. An interesting 
and plausible possibility comes from the fact that the gluon propagator may be non- 
perturbative and IR singular because of the confining nature. If this is the case, such 
IR singular interaction would induce an inhomogeneous chiral condensate as elucidated 
in section 7.1 or the quarkyonic chiral spiral [214]. 

8.2. QCD ai Ac = 2 

A QCD-like theory with Ac = 2 [245, 246] is also an interesting limit opposite to the 
case of quarkyonic matter; small Ac instead of large Ac. Two-colour QCD is free from 
the sign problem if the number of degenerate flavours is even, so that Monte-Carlo 
simulation at finite baryon density is possible [247]. Two-colour QCD shares many 
non-perturbative features with Ac = 3 QCD such as confinement, chiral symmetry 
breaking and superfluidity. 

In two-colour QCD, baryons consist of quark pairs qq which are bosons. Therefore, 
nuclear matter composed of fermionic baryons in Ac = 3 is replaced by a superfluid 
state of bosonic baryons in Ac = 2. The order parameter of superfluidity is the 
colour-singlet baryon condensate (qq) . To investigate the phase structure of two-colour 
QCD, such baryon condensate has been computed both in numerical simulations and 
in analytical strong-coupling expansion [248, 249, 250]. There is also a suggestive 
result that supports the idea of quarkyonic matter in Ac = 2 from lattice simulations 
[251], which is consistent with the model analysis [252]. The correlation functions 
and excitation spectra can be studied both analytically and numerically for Ac = 2. 
In particular, the low-energy chiral Lagrangian at finite baryon density as well as 
Leutwyler-Smilga type spectral sum rules have been derived [253] . In- medium hadron 
spectra as a function of the baryon chemical potential have been also investigated 
[254, 255, 256]. As can be seen in these examples, two-colour QCD continues to be a 
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Table 4. Correspondence between the boson-fcrniion mixture in ultracold atoms 
(such as '^'^Rb and '^''K mixture) and the diquark-quark mixture in 2-flavour QCD. 



valuable testing ground for studying hot/dense QCD at Nc — 3. 

8.S. Ultracold atoms 

High density QCD matter and ultracold atomic systems, although differing by some 
twenty orders of magnitude in energy scales, share analogous physical aspects [257] . 
Phenomenological studies of QCD indicate a strong spin-singlet diquark correlation 
inside the nucleon [258]. Therefore, it may be a good starting point to model the 
transition from 2-flavour quark matter at high density to nuclear matter at low density 
in terms of a boson-fermion mixture, in which small-size diquarks are the bosons, 
unpaired quarks the fermions and the extended nucleons are regarded as composite 
Bose- Fermi particles [259]. 

Recent advances in atomic physics have made it possible indeed to realize a boson- 
fermion mixture in the laboratory. In particular, tuning the atomic interaction via a 
Feshbach resonance allows formation of heteronuclear molecules, as recently observed 
in a mixture of ^^Rb and ^°K atomic vapours in an optical dipole trap [260]. An 
analysis of such non-relativistic mixture of cold atoms indicate that the BCS-like 
superfluidity of composite fermion (N=b-|-f ) with a small gap is a natural consequence 
of the strong boson-fermion attraction [259], which may explain the reason why the 
fermion gap in nucleon superfluidity can be an order of magnitude smaller than the gap 
in colour superconductivity. A possible correspondence between cold atoms and QCD 
is summarized in table 4. Fuller understanding, both theoretical and experimental, 
of the boson-fermion mixture as well as a mixture of three species of atomic fermions 
[261, 262] may further reveal properties of high-density QCD. 



9. Summary and concluding remarks 

In this article we reviewed the current status of theoretical investigations to explore the 
QCD phase structure at finite temperature (T) and finite baryon chemical potential 
(hb)- There are (at least) three fundamental states of matter in QCD: the hadronic 
matter with broken chiral symmetry and quark confinement in the low-T and low-/iB 
region, the quark-gluon plasma at high T and the colour superconductivity at low T 
and high fiB- On top of these states, some exotic phenomena have been conjectured, 
e.g. the chiral-density waves, the crystalline colour superconductivity, the gluonic 
phase, the quakyonic matter and the quark-hadron continuity, as covered in this 
article, and even more possibilities are still developing. Lattice QCD simulations, the 
Ginzburg-Landau- Wilson approach and effective theories of QCD are useful theoretical 
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tools to study these phenomena. 

From the experimental point of view, the QCD phase transitions at high T with 
Hb/T < 1 can be studies by using high-energy heavy-ion collisions at RHIC and 
LHC. The QCD phase diagram at relatively low T with hb/T >^ 1 may also be 
probed in the future facilities with lower-energy heavy-ion beams. Besides, recent 
attempts to determine the mass and the radius of neutron stars from X-ray bursts 
would be an alternative way to access the equation of state of dense QCD [263, 264]. 
Furthermore, gravitational waves [265] and neutrinos from supernova explosions [266] 
or cooling processes in the neutron star [267] carry valuable information of dense QCD 
matter. Dense QCD shall continue to be one of the most fascinating theoretical and 
experimental topics in particle, nuclear and astro physics. 
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